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We report multi-component high-entropy materials as extraordinarily robust diffusion barriers and clarify 
the highly suppressed interdiffusion kinetics in the multi-component materials from structural and 
thermodynamic perspectives. The failures of six alloy barriers with different numbers of elements, from 
unitary Ti to senary TiTaCrZrAlRu, against the interdiffusion of Cu and Si were characterized, and 
experimental results indicated that, with more elements incorporated, the failure temperature of the 
barriers increased from 550 to 900 °C. The activation energy of Cu diffusion through the alloy barriers was 
determined to increase from 1 10 to 163 kj/mole. Mechanistic analyses suggest that, structurally, severe 
lattice distortion strains and a high packing density caused by different atom sizes, and, thermodynamically, 
a strengthened cohesion provide a total increase of 55 kJ/mole in the activation energy of substitutional Cu 
diffusion, and are believed to be the dominant factors of suppressed interdiffusion kinetics through the 
multi- component barrier materials. 

In 2004, Yeh reported a brand-new concept of alloy design, named high-entropy alloys (HEAs), with equimolar 
incorporations of multi-principal (more than five) metallic elements\ Thermodynamically, owing to high 
mixing entropies, simple solid-solution structures rather than any complex phases or intermetallics are 
typically formed in HEAs^'^; structurally and kinetically, severe lattice distortions caused by different atom sizes 
suppress diffusion and nucleation/growth rates, generally yielding nanocrystalline or even amorphous struc- 
tures'*. In the past years, HEAs have been intensively studied, and many extraordinary properties including a good 
mechanical performance^, a unique elasticity^, high plasticity^, and high saturation magnetization and electrical 
resistivity^ been discovered. Films of HEA nitrides (HEANs), such as (AlCrTaTiZr)N, with simple solid-solution 
or nanocomposite structures and very good mechanical properties (high hardness and wear resistance) were 
afterwards developed as protective hard coatings^'*^. Attributed to their high thermal stability, thin HEA and 
HE AN layers have also been considered for further applications to diffusion barriers *\ Robust diffusion barriers 
are of importance for use as interfacial layers particularly in the interconnects of integrated circuits (ICs), e.g. 
between Cu wires and dielectrics or between metal pads and solders, to inhibit the interdiffusion of adjacent 
materials and the formation of compounds (e.g. silicides)*^. Refractory transition metals, such as Ti and Ta*^, and 
the nitrides of unitary transition metals, including TiN and TaN of strong ionic bonds *^'*^, were first developed 
and practically used in interconnects, but however unsuited for applications to next- generation interconnects due 
to their high solid solubility to Cu and their columnar structures with boundaries as fast diffusion paths. Thinner 
and more effective barriers have been investigated in recent years, which typically comprise: (1) ternary compo- 
nents with an amorphous structure to diminish diffusion paths, such as Ru-Ti-N and Ru-Ta-N*^"*^, or (2) layered 
structures with interface mismatches to elongate diffusion paths, such as Ru/TaN and Ru/TaCN*^"^°. 
Nevertheless, in the IC manufacturing generation below 20 nm, robust and ultrathin barrier layers (<3 nm) 
with a high diffusion resistance are strongly demanded but rarely reported. Though self-forming (barrier-less 
alloy) and self- assembled (organic monolayer) diffusion barriers have been developed in the past few years^*'^^, 
their thermal stability is yet concerned. 

Recent studies indicated that, by the multi- component-induced lattice distortions and the layered structure- 
caused interface mismatches, HEA and HE AN barriers such as senary AlCrRuTaTiZr and (AlCrRuTaTiZr)N and 
their stacking structures with a thickness down to only 4 nm presented an extreme resistance to the interdiffusion 
of Cu and Si^^~^^. Their high endurance temperature (almost 900°C) against interdiffusion suggests the high 
potential of HEAs and HEANs as very promising barriers for practical applications to Cu interconnects. 
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Unfortunately, except a general concept of "lattice distortions", not 
an exact mechanism why HEAs and HEANs provide a high resist- 
ance to interdiffusion has been established. An amorphous structure 
is not simply the answer because thin sputtered Ti and TiN films are 
also possibly amorphous but do not present such a high resistance. 
The question remains. Hence in this study, six alloy films with dif- 
ferent numbers of metallic elements (nA), from unitary Ti (lA), 
binary TiTa (2A), ternary TiTaCr (3A), quaternary TiTaCrZr 
(4A), quinary TiTaCrZr Al (5 A) to senary TiTaCrZr AlRu (6 A) (ele- 
ment added following the development sequence of diffusion bar- 
riers), were prepared and characterized; the resistance of ultrathin 
alloy diffusion barrier layers (DB^aj ~7 nm thick) to the interdiffu- 
sion of Cu and Si (in term of the failures of Si/DB„A/Cu film stacks), 
along with the interdiffusion kinetics through the alloy barriers, were 
determined to examine the enhanced resistant ability with more 
metallic elements incorporated. Mechanistic studies of structural 
and thermodynamic factors, from the aspects of lattice distortion, 
cohesive energy and packing density, were also performed to clarify 
the suppressed interdiffusion kinetics in multi-component high- 
entropy alloys. 

Results 

Compositions and structures. Compositional and structural 
analyses, as shown in Figure 1, indicate the uniform and simple 
solid- solution structures of the six alloy films investigated in the 
present study. From the elemental contents of the as-deposited 
thick alloy (nA) films given in Figure 1 (a), it was known that the 
metallic elements added in the alloys were in near-equimolar ratios, 
only except that in the binary alloy the Ta content was unexpectedly 
higher than Ti (Ta : Ti ~ 3 : 1) due to the resputtering of the light Ti 
compared to the heavy Ta. A few content (~3 at.%) of oxygen 
existed owing to the limited vacuum of the sputtering chamber. 
Although the elemental contents of the thick alloy films do not 
certainly equal to the compositions of thin diffusion barrier layers 
because they may slightly vary with film thickness, the measured data 
are still believed to represent the good composition uniformity of the 
thin layers. The XRD pattern of each thick alloy film in Figure 1 (b) 
presented only one set of peaks, revealing the simple solid-solution 
structures of the alloys even with the additions of multiple 
components. For the unitary alloy, main diffraction peaks at 20 = 
35.0°, 38.3° and 40.1° (corresponding to the (100), (002) and (101) 
lattice planes of Ti; JCPDS #44-1294) indicated a hexagonal close- 
packed (hep) crystal structure, while for the binary alloy, the peaks at 
38.4°, 55.2° and 69.5° (to (110), (200) and (211) planes) differently 
suggested a body-centered cubic (bcc) structure because a large 
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Figure 1 | (a) Elemental contents and (b) XRD patterns of as-deposited 
packed, bcc: body-centered cubic). 



portion of bcc Ta (JCPDS #04-0788) was added. In comparison for 
the ternary to senary alloys, broad hills, rather than sharp peaks, were 
observed, indicating the appearance of amorphous structure. 
Nevertheless, the positions (2^) of the hill tops around 38.5-40.5° 
and 69.8-70.9° (sin^ ratio ~ 1 : V^) well accorded with the (110) 
and (211) peaks of bcc structure, suggesting a tendency of forming 
a bcc-based, short-range-ordered structure (embedded in the 
amorphous matrix) in these alloys. Peak shifts were noticed, which 
accorded with the variations in the average atom sizes (following the 
Bragg law) of the alloys with more metallic elements incorporated, as 
Usted in Supplemental Table SI. 

Diffusion resistance. The resistance of thin barrier layers (DB„a) 
with different numbers of metallic elements to the interdiffusion of 
Cu and Si is determined from the failures of the SI/DB^a/Cu film 
stacks, which are characterized from SEM morphologies, XRD 
patterns, electrical resistivities, and also TEM microstructures 
before and after thermal annealing, as presented in Figures 2 to 6. 
Figure 2 first shows the SEM morphologies of the as-deposited and 
thermally annealed SI/DB^a/Cu stacks. As the annealing 
temperature increased to the endurance limits where the diffusion 
barriers still resist obvious interdiffusion of Cu and Si, the grains of 
the top Cu films grew, accompanied by slight Cu agglomerations to 
form voids or even islands (in particular for the SI/DBia/Cu stack at 
500°C and SI/DBsa/Cu at 600°C), but still exhibited smooth surfaces. 
The stacks failed at extreme temperatures (defined as the failure 
temperatures) once Cu penetrated through the barrier layers into 
and reacted with Si to form typical particulate-like Cu silicides. 
With more metallic elements incorporated in the alloy barriers, 
clearly, the failure temperature increased from 600 to 900 °C. 

The XRD patterns and electrical resistivities of the as-deposited 
and thermally annealed SI/DB^a/Cu stacks, as plotted in Figures 3 (a) 
and (b), respectively, showed the same failure trend with the number 
of metallic elements. In the XRD pattern of the as-deposited stacks, a 
slightly broad peak at 43.3° corresponded to the (111) lattice plane of 
face-centered cubic (fee) Cu with a fine grain size. For the annealed 
stacks at the endurance limits, the Cu (1 1 1) peaks sharpened because 
of grain growth, while at their failure temperatures, two additional 
peaks at 44.7° and 45.3° corresponding to the (320) and (312) lattice 
planes of typical CusSi compound appeared, revealing the severe 
interdiffusion of Cu and Si. The formation of Cu silicides led to 
consequent drastic increases in the electrical resistivities of the failed 
stacks. Originally, high electrical resistivities of the as-deposited 
stacks were measured, above 10 |iQ-cm, due to the defects in the 
Cu films. With thermal annealing till the endurance limits, the 
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Figure 2 | SEM surface morphologies of as-deposited and thermally annealed SI/DB^a/Cu film stacks {n: number of metallic elements; Cu film 
agglomerations into islands on the film stacks with DBia at 500°C and DB2A at 600°C); insets: (up left) as-deposited with a typical Cu film surface, (down 
right) failed with a typical silicide formation. 



electrical resistivities gradually decreased to only 2.5 |iQ-cm, close to 
the value of pure Cu, owing to Cu grain growth and defect elimina- 
tions. Then at the failure temperatures (600 to 900 °C for SI/DBia/Cu 
to Si/DB6A/Cu respectively), a jump in resistivity to a high value 
occurred as a consequence of the abundant formation of Cu silicides 
(note: the resistivities of some regions on Si/DBjA/Cu and Si/DB2A/ 
Cu stacks early jumped at 500 and 600°C, respectively, owing to Cu 
agglomerations to discontinuous islands). In comparison, Si/Cu 
without a barrier failed very early at 400-500°C similar to the 
reported temperature^^. 

The above analyses verified that the SI/DBia/Cu film stack com- 
pletely failed at 600°C, while Si/DBgA/Cu did not until 900°C. 
However, it was noted that, by careful examinations of magnified 
SEM morphologies, XRD patterns and electrical resistivities 



(Supplemental Figure SI), the SI/DBia/Cu and SI/DBsa/Cu stacks 
might have possibly failed at 500-600°C and 800-900°C, respect- 
ively, because very small compound particles formed, very low (320) 
and (312) peaks of sihcides appeared, and electrical resistivities 
slightly rose. As plotted in Figure 4, it is concluded that the failure 
temperature of the film stacks, namely the diffusion resistance of the 
barriers, increased with the number of incorporated metallic ele- 
ments, from 500-600°C (noted as 550°C) for the unitary, 700°C 
for the binary, 800°C for the ternary and quaternary, 800-900°C 
(850°C) for the quinary, to 900°C for the senary alloy. 

The TEM cross-sectional microstructures and lattice images of the 
as-deposited and the 400-600°C annealed Si/DBiA/Cu and 700- 
900 °C annealed Si/DB6A/Cu film stacks shown in Figure 5 confirmed 
the markedly enhanced diffusion resistance of multi-component 
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Figure 3 | (a) XRD patterns and (b) electrical resistivities of as-deposited and thermally annealed SI/DB^a/Cu film stacks ( n: number of metallic elements; 
AD: as-deposited; OA: Si/Cu without a diffusion barrier; following "lA" and "2A": Cu film agglomerations). 
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Figure 4 | Failure temperatures of Si/DB^A/Cu film stacks (n: number of 
metallic elements). 

alloys. In the as-deposited stacks, both the DBia and DB6a layers of 
only 7 nm thick showed an amorphous structure, while the Si sub- 
strate and the Cu film (interplanar spacings 0.209 and 0.181 nm, 
(111) and (200) lattice planes; JCPDS # 85-1326) exhibited a crystal- 
line, ordered structure. Between each two layers, a clear, continuous 
interface was identified. Compared to the thick unitary alloy film of a 
crystalline hep structure and the thick senary alloy film of a bcc- 
based, short-range- ordered structure, the as-deposited ultrathin bar- 
rier layers more easily formed an amorphous structure^^'^^. However, 
a possibility should still be noted that a short-range ordering might 
exist in the amorphous (disordered) matrix^^, which needs to be 
further investigated. After annealing of SI/DBia/Cu at 400 to 
500°C and Si/DBgA/Cu at 700 to 800°C, clear interfaces and the 
amorphous structure of the alloy barriers basically remained despite 
few indistinct short-range orderings (several atomic layers in size) 
that might have possibly formed in the amorphous matrix of the 
500°C-annealed DBia and the 800°C-annealed DBgA layers. 
However, in the 500°C-annealed SI/DBia/Cu, few CuTi3 nanocrys- 
tallites (interplanar spacing 0.227 nm, (114) plane; JCPDS #50-1476) 
were noticed, suggesting the interdiffusion of Cu and Ti as well as the 
beginning of failure. The SI/DBia/Cu completely failed at 600°C, and 
the Si/DBgA/Cu did at 900°C; diffused interfaces and CuTi3, Cui5Si4 
(interplanar spacing 0.340 nm, (220) plane; JCPDS #76-1800) and 



CusSi (interplanar spacing 0.280 nm, (210) plane; JCPDS #04-0841) 
nanocrystallites were found. 



Activation energy of Cu diffusion. From Pick's second law^°, for 
non-static diffusion in a semi-infmite system which is similar to the 
interdiffusion of Cu and Si in Si/DB„A/Cu film stacks through the 
thin alloy barrier layers, the concentration C of diffusing species as a 

M 

function of distance x and time t, is given as C = —= exp 



/nDt ^ V 

where M and D are the thickness and diffusivity of the diffusing 
species, respectively, and the diffusion distance (at a constant 
concentration) is given as x = 2\/Dt. Accordingly, to determine the 
diffusion distances and diffusivities of Cu through the alloy barriers, 
the concentration distributions of Cu and Si elements in the as- 
deposited and thermally annealed SI/DBia/Cu and Si/DB6A/Cu 
stacks at different temperatures were detected by TEM line 
scanning, as typically presented in Supplemental Figure S2. The 
curve fittings of the as-received data were performed using a 
sigmoid function to yield smooth elemental concentration profiles, 
as plotted in Figure 6 (the interfaces between the barriers and the 
native oxide were offset to zero in distance). Obviously, as 
temperature increased, Cu penetrated far into Si substrates 
through the alloy barriers. As presented in Figure 7 (a), the 
diffusion distances of Cu, at a fixed concentration of 50 at.%, in 
the annealed SI/DBia/Cu and Si/DB6A/Cu stacks (referred to the 
as-deposited stacks) normally increased with temperature but 
decreased with the number of incorporated metallic elements. The 
diffusivities of Cu at different temperatures were further determined 
by using x = 2\/Dt, as plotted in Figure 7 (b) (vs. 1 /temperature). In 
the SI/DBia/Cu stack the diffusivity of Cu reached 4.1 X 10"'° mVs 
at 600°C, while in the Si/DBgA/Cu it was only 1.6 X 10""^ mVs at 
700°C and 8.0 X 10"'° mVs at 900°C, obviously indicating the 
suppressed interdiffusion by the multi- component barriers. 

From the experimentally measured logarithmic diffusivity- to- 1/ 
temperature relations in Figure 7 (b) and the equation 

D = Do exp ^ — in which the diffusivity constant Dq is a fre- 
quency factor, R the gas constant, and Tthe temperature^", the activa- 
tion energy of Cu diffusion in the SI/DBja/Cu and Si/DBeA/Cu 
stacks, through the unitary and senary alloys, is then determined 
as 1 10 and 163 kj/mol, respectively, as given in Figure 8. The increase 
in activation energy, 53 kJ/mol, yields a lowered diffusivity of Cu and 
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elucidates the markedly enhanced diffusion resistance of the senary 
alloys. Further according to the measured failure temperatures of Si/ 
DB„a/Cu stacks plotted in Figure 4, the activation energy of Cu 
diffusion through the binary to quinary alloys is proportionally inter- 
polated, as also given in Figure 8 that clearly reflects the raised dif- 
ficulty of Cu diffusion through the alloy barriers with more metallic 
elements incorporated. 

Discussion 

Several models for Cu diffusion in conventional barrier materials, 
established by using molecule dynamic simulations, have been 
reported^^"^^, including the diffusion in metals via vacancies along 
different directions and that in nitrides via vacancies or interstitial 
sites. The model for diffusion in bcc Ta via vacancies suggests the 
smallest activation energy of diffusion along the [111] direction 
(0.65 eV, 63 kj/mol), but higher energy along [100] (2.7 eV, 
261 kJ/mol) and [110] (5.8 eV, 560 kJ/mol)'\ In the present study, 
the similar-level but ascending activation energy with increasing the 
number of metallic elements is believed to reveal that Cu diffusion in 
the alloy barriers proceeds substitutionally via vacancies along the 
most facile direction but is possibly hindered by some factors dis- 
cussed below. For the substitutional diffusion of atom from a lattice 
site to another, two requirements need to be fulfilled: (1) an adjacent 
vacancy (activation energy of vacancy formation, AGy) and (2) a 
bond breaking and an atom transfer (activation energy of movement. 



AG^), both of which compose the activation energy of diffusion^°. 
Accordingly, three important factors for the increased activation 
energy of Cu diffusion, i.e. suppressed interdiffusion of Cu and Si, 
in the present multi- component alloy barriers are considered to be, 
structurally, the strain energy of lattice distortions and the packing 
density of atoms, and, thermodynamically, the cohesive energy of 
atoms, all of which concern vacancy formation and atom movement. 

(1) Lattice distortion strain energy. Figure 9 (a) schematically 
illustrates the normal (regular, isotropic) lattice structure of 
unitary alloy and, differently, the distorted (atom displaced and/or 
bond shortened/stretched) structure of multi-component soUd- 
solution alloy (random distribution with different atom sizes and 
cohesion with different bond strengths). In the unitary lattice with 
the same stress state in all directions, atoms diffuse via adjacent 
vacancies along the most facile direction (e.g. [Ill] in a bcc 
crystal) of the lowest energy barrier^\ However, in the multi- 
component lattice with different stress distributions in different 
directions, the lattice strain energy is anticipated to influence the 
site of vacancy formation and the direction of atom movement; 
vacancies will form and atoms will migrate more frequently along 
more distorted directions (not certainly [111]). The total diffusion 
paths will be elongated, and the average difficulty in diffusion (i.e. the 
activation energy of diffusion) will increase, lowering the diffusivity 
of atoms in the multi -component alloys. The lattice distortions 3 in 
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the multi- component alloys can be calculated by using the equation 

d = ^y^X/ [\ — 3^ in which X/ is the fraction of the /th 

component, and and r the /th and the average atom radii, 
respectively^^. The lattice distortion strain energy Ku strain is then 
theoretically calculated from the equation /SiU strain = {^^){^^) d 
where E is the average modulus (listed in Supplemental Table SI) 
and d the average displacement of atoms {d = 3f), as plotted in 
Figure 9 (b) with the number of incorporated metallic elements. 
Compared to the zero energy of the unitary alloy without distor- 
tions, the lattice strain energy of the present multi-component 
alloys increases to 16.8 kj/mol with four metallic elements and 
slightly decreases to 11.9 and 14.6 kJ/mol with five and six 
elements, respectively. The change in the electrical resistivities of 
thick alloy films, as given in Supplemental Figure S3, confirms the 
variation of lattice distortions with the number of metallic elements. 
Besides a typical film defect-caused increase in experimentally 
measured electrical resistivity from a theoretically calculated value 
(by the rule of mixture) an extra increase (delta, relative to the 
typical increase of unitary alloy without lattice distortions) induced 
by lattice strains follows the variation trend of lattice distortions 
though a slight inconsistency exists for the quaternary and quinary 
alloys. The high lattice strain energy in the multi- component alloys 
with different- size atoms is also a dominant factor for the high 
stability of their amorphous structure in intense annealing, as sug- 
gested in Results that the temperature for the amorphous senary alloy 
to endure crystallization, 800 °C, is much higher than that for the 
unitary alloy, only 500°C. Crystallization of multi-component alloys 
will cause much severer lattice distortions and higher lattice strain 
energy, while a topological instability will facilitate an amorphous 
solid solution, especially at high temperatures where an entropic 
force assists the formation and stabilization of solid-solution 
phases (random distributions of different- size atoms) 

(2) Cohesive energy. As also seen in Figure 9 (a), in the unitary lattice 
with the same cohesion between any adjacent atoms, atoms move via 
adjacent vacancies along the most facile direction^ \ In the multi- 
component lattice with anisotropic cohesions in different direc- 
tions, however, high cohesive energy is expected to inhibit a bond 
breaking for vacancy formation and atom movement, and to change 
the direction of atom movement; vacancies will form and atoms will 
migrate more frequently along less strengthened directions. The 
average bond breaking energy will increase, the total diffusion 
paths will be elongated, and the average difficulty in diffusion (i.e. 
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the activation energy of diffusion) will increase, then lowering the 
diffusivity of atoms in the multi-component alloys as well. Supple- 
mental Table S2 lists the mixing enthalpies AH^^^ and cohesive 
energy Qy of any two elements in the present multi- component 
alloys where the cohesive energy of the /th and jth components is 

obtained by Miedema's thermodynamic model Qij = " ^ — - — 



^^m/x 37-39 



. From Takeuchi's regular- solution model Q!^^^ = 



QijXiXf^, the cohesive energy Q"'''' of the present multi- 

j=i+i 

component alloys is further theoretically calculated. Relative to the 
cohesive energy of the unitary alloy, 14.2 kJ/mol, the average 
increases in the cohesive energy, AH'"''', of the multi-component 
alloys are found to increase to 38.1 kJ/mol with the number of 
metalUc elements, as plotted in Figure 9 (b). 

(3) Packing density. The packing density of atoms in a unitary alloy 
with a bcc or fcc/hcp structure is 0.68 or 0.74, respectively. How- 
ever, in a multi- component solid- solution alloy that is mixed with 
different- size atoms, the packing density will increase^\ The shrink- 
age in free volume is believed to diminish vacancy formation/con- 
centration, additionally increasing the difficulty in diffusion (i.e. the 
activation energy of diffusion) and lowering the diffusivity of atoms 
in the multi-component alloys. From Danisch's random-packing 
model of different- size balls^\ the packing densities of multi- 
component alloys with different coordination numbers and atom- 
r 

size ratios (z 1 + ^, ^: lattice distortions) are theoretically calculated, 

as plotted in Supplemental Figure S4. For the hep-based unitary alloy 
(with an atom-size ratio of 1), the packing density S is normally 0.740, 
while for the bcc-based binary (an atom-size ratio of 1.006) to senary 
(a ratio of 1.074) alloys, the size-dependent densities S' increases 
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(from an original S of 0.680) to 0.683 and to 0.711, respectively. 
Combining the contributions of the lattice distortion strain energy, 
the increase in cohesive energy and the multiplied packing density, 
the total increases in the activation energy of diffusion in the alloys 

are predicted as AQ(i) = [AUstrain + AH"^^^) — , as plotted in Figure 9 

(b). For the senary alloy, the calculated energy increase, 55 kj/mol, 
nearly equals the measured increase, 53 kJ/mol (from 110 to 163 kj/ 
mol), suggesting the validity of the prediction. 

Figure 10 comparatively presents the experimentally measured 
activation energy of Cu diffusion and theoretically calculated values 
(110 kJ/mol offset, referred to the experimental value of the unitary 
alloy) in the present alloys with different numbers of metallic ele- 
ments. Similar trends of the experimental and calculated data with 
the number of metallic elements are observed though, small differ- 
ences between them are found yet, possibly attributable to a 
'Volume-to-area" effect. The reported models for Cu diffusion in 
alloys basically consider atom motion on a two-dimensional 
^j.g^3i-33. however, the present calculations are on the basis of 
three-dimensional volume. By considering a volume-to-area trans- 
formation, the theoretical calculations are thus modified as 

AQ(2) = /:(AQ(i))^^^ {k =^ 3.7), as also plotted in Figure 10. The 
modified calculations well match the experimental values, conclud- 
ing the two-dimensional effects of the dominant three-dimensional 
structural (distortion/packing) and thermodynamic (cohesion) fac- 
tors on the suppressed interdiffusion kinetics in the multi- compon- 
ent high-entropy alloys. 

In summary, six solid-solution alloy barriers, from unitary Ti to 
senary TiTaCrZrAlRu, with different numbers of metallic elements 
of different atom sizes and cohesive energy were characterized. As 
experimentally determined, with more metallic elements added, the 
failure temperature of the alloy barriers against the interdiffusion of 
Cu and Si increased from 550 to 900 °C, and the activation energy of 
Cu diffusion in the alloys increased from 110 to 163 kJ/mole, indi- 
cating an improved diffusion resistance of the multi-component 
high-entropy alloys. Mechanistically, the substitutional diffusion of 
Cu atoms via vacancies, with the consideration of a volume-to-area 
transformation, is expected. Dominate factors that provide a theor- 
etical increase of 55 kJ/mol in the activation energy of Cu diffusion in 
the present multi-component alloys and suppress the interdiffusion 
kinetics of Cu and Si include, structurally, a high lattice distortion 
strain energy and a high packing density, and, thermodynamically, a 
high cohesive energy, to diminish the formation of vacancies and to 
hinder the movements of atoms, yielding an extreme diffusion 
resistance. 



Methods 

Thin and thick alloy films with different numbers of metallic elements (nA), including 
unitary Ti (lA), binary TiTa (2A), ternary TiTaCr (3A), quaternary TiTaCrZr (4A), 
quinary TiTaCrZrAl (5A) and senary TiTaCrZrAlRu (6A), were deposited by typical 
magnetron sputtering and investigated. Thin alloy diffusion barrier layers (DB„a> 
~7 nm) were prepared for examining their resistance to the interdiffusion of Cu and 
Si; thick alloy films {nA, 600-1000 nm) were for characterizing compositions, crystal 
structures and electrical resistivities. Sputtering targets were fabricated by arc-melting 
equimolar constituent elements in vacuum, and then cutting and polishing the 
solidified bulks into discs with a diameter of 50 mm. The thin barrier layers and thick 
alloy films were deposited on Si substrates at a radio -frequency power of 150 W and a 
substrate bias of — 100 V, at room temperature, in an Ar atmosphere (gas flow rate of 
30 seem, working pressure of 6 X 10"^ torr). A top Cu film (50 nm) was deposited on 
the thin barrier layers at a power of 50 W (Ar flow rate of 13 seem, working pressure 
of 5 X 10"^ torr) to form Si/DB„A/Cu film stacks. Thermal annealing of the film stacks 
was then carried out in a tube furnace (base vacuum of 4 X 10"^ torr) at 300 to 900° C 
(ramping rate of 5°C/min) for 30 minutes in an N2/H2 atmosphere (total gas flow rate 
of 500 seem, H2 flow ratio of 5%; reducing atmosphere to prevent the film stacks from 
oxidation) to investigate the resistance of the barrier layers to the interdiffusion of Cu 
and Si. For comparison, Si/Cu without applying a diffusion barrier layer was also 
studied. 

The elemental compositions of the thick alloy films were determined by electron 
probe micro -analyses (EPMA, JEOL JXA-8800M). A scanning electron microscope 
(SEM, JEOL JSM-6700F) was used to observe the surface morphologies of as- 
deposited and thermally annealed Si/DB„A/Cu film stacks, and a glancing incident 
angle (0.5°) X-ray diffractometer (XRD, Rigaku Dmax 2000) was applied to analyze 
the crystal structures of the thick alloy films and the as -deposited and thermally 
annealed Si/DB„A/Cu stacks. The sheet resistance of the thick alloy films and the film 
stacks was measured by using a four-point probe station (Keithley 2400) to calculate 
electrical resistivities. The microstructures and lattice images of the as-deposited and 
thermally annealed Si/DB„A/Cu film stacks (with DBia and DBsa) were observed by a 
high-resolution transmission electron microscope (TEM, JEOL JEM-2100F), and the 
elemental concentration profiles of Cu and Si in the film stacks were determined by 
line scanning analyses using energy dispersive spectrometry (EDS, Oxford). 
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